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Abstract 

We generalise the description of the dynamics of the order book of financial markets in terms of a 
Brownian particle embedded in a fluid of incoming, exiting and annihilating particles by presenting 
a model of the velocity on each side (buy and sell) independently. The improved model builds on the 
time-averaged number of particles in the inner layer and its change per unit time, where the inner 
layer is revealed by the correlations between price velocity and change in the number of particles 
(limit orders). This allows us to introduce the Knudsen number of the financial Brownian particle 
motion and its asymmetric version (on the buy and sell sides). Not being considered previously, the 
asymmetric Knudsen numbers are crucial in finance in order to detect asymmetric price changes. 
The Knudsen numbers allows us to characterise the conditions for the market dynamics to be 
correctly described by a continuous stochastic process. Not questioned until now for large liquid 
markets such as the USD/JPY and EUR/USD exchange rates, we show that there are regimes 
when the Knudsen numbers are so high that discrete particle effects dominate, such as during 
market stresses and crashes. We document the presence of imbalances of particles depletion rates 
on the buy and sell sides that are associated with high Knudsen numbers and violent directional 
price changes. This indicator can detect the direction of the price motion at the early stage while 
the usual volatility risk measure is blind to the price direction. 
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I. INTRODUCTION 


The mean free path length, the average distance travelled by a moving particle between 


successive collisions, is an important quantity to characterise fluids [lj. In statistical physics, 
especially hydrodynamics, the Knudsen number ( 2 - 4], which is defined as the ratio of the 
mean free path length to a representative particle length scale, is used to determine whether 
the system can be described in the continuum limit or needs a description accounting for 
discrete or particle effects. Generally, the fluid can be well approximated by continuous 
mathematical equations when the Knudsen number is below 0.1, so that sufficiently many 
collisions occur within a particle size to erase the influence of any discreteness. Thus, 
for a Knudsen number below 0.1, one uses in general the Navier-Stokes equation, while 
the Boltzmann equation is needed for Knudsen number larger than 1, with a complicated 
transition in between. There are long standing problems characterised by high Knudsen 
numbers, including dust particle motions through the lower atmosphere, satellite motion 
through the exosphere, microfluidics js| and micromechanical systems. These “high Knudsen 
number flows” possess novel properties not shared by fluid flows at low Knudsen numbers, 
which are being investigated with novel instruments and theoretical methods jfj-fl. 

Beyond the scale of the mean free path, particles follow random paths, a large class of 
which are called Brownian motions. This later concept describes the random motion of small 
objects in fluid media, which was discovered by Jan Ingenhousz in 1784 when observing the 
irregular motion of coal dust particles on the surface of alcohol 10] and Botanist R. Brown in 
1827 when observing pollen grains under a microscope. The mechanism of this phenomenon 
was finally clarified by independently Einstein and Smoluchowski, explaining the random 
motion of the Brownian particle as due to the incessant collisions with the embedding fluid 
particles, proving theoretically the particle nature of matter jll, 12]. This idea is one of the 
fundamental starting points of statistical physics, which connects micro-scale to macro-scale 
phenomena. 

Several years before Einstein’s random walk theory, Bachclier introduced the random 
walk model to describe the dynamics of prices of assets traded in stock markets [3]. Indeed, 
prices of traded assets fluctuate incessantly and randomly up and down as a result of the 
aggregation of all traders’ decisions. The concept of random walks to characterise the 
aggregate actions of motivated and diligent investors was a significant breakthrough that 
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has further developed into the “efficient market hypothesis” (see M for a short history and 
review), which forms the foundation of financial engineering [lfl. Finance and physics are 
both founded on the theory of random walks and their many generalisations, going back 
to more than one hundred years of intertwined history [14J]. However, the meaning of the 
Knudsen number in finance has not be elucidated until now and this will be the principal 
goal of this article. 

In 1963, Mandelbrot documented that the scale of volatility (defined as the absolute price 


changes in a given time interval) is nonstationary and t 
change is asymptotically described by a power law tail 


re probability distribution of price 


161] . These two properties are now 


recognised as fundamental stylised facts of financial price series. Building on analogies with 
the many power law distributions in physics, this has later excited the interest of statistical 
physicists to analyze financial time series further, using the large amount of available data. 
Recently, several novel properties of financial prices series have been documented, such 
as the negative auto-correlation of price changes at short times, the abnormal diffusion 
of prices BA, the long-auto correlation of volatility j^ioj], the long-memory process of 


sign of orders 


21 


22], the property that implicitly underlies the present considerations as 


23 


24], large price changes characterized by the gap of a limit 


financial multifractality 
order book (containing no quotes between prices) 25|, the existence of endogenous feedback 
mechanisms that are well characterised by the self-excited Hawkes model 26[, rich market 


impact dynamics revealed in nonlinear price changes caused by submitted orders 
and the Fokker-Plank description for the queue dynamics of the order book [fiol ]. 

In this paper, we analyze financial markets from the micro scale view of traders 


29j, 


33|], 


We focus on a foreign exchange market, which is the place where a trader can exchange 
his currency to another currency via a trade with another trader. As an example, in the 
USD-JPY currency market, we call a trader who wants to change JPY to USD a buyer and 
USD to JPY a seller. There are two possible actions available to each trader. The first one 
is to submit a market order that is immediately executed at the best price in the market 
at the time of the submission. The second possible action is a limit order that enables a 
trader to book his order in the market with a specific price and quantity, and which can 
be canceled later if he wishes. If two limit orders have been placed at the same price, the 
tie is broken on a first-come first-served basis. The aggregation of all limit orders is called 
the order book, in which the limit orders both on the buyer side (bids) and on the seller 
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side (asks) are placed on the price axis while the volume associated to each limit price is 
represented as a 2-D graph to be shown later. A trade occurs when either a market order or 
a limit order meets the demand/supply of existing limit orders. The prices and volumes in 
the order book and thus the realised market price time series result from the aggregation of 
all traders’ decision making processes. Describing in details the full dynamics of the order 
book process is a challenging task that can benefit from the perspective of statistical physics 
applied to financial markets. 

Recently, the celebrated fluctuation-dissipation Theorem (FDT), which states that the 
response of a system in equilibrium to a small applied force is the same as its response to a 
spontaneous fluctuation [3J,|35[, has been found to hold for the average dynamics of the limit 
order prices in the order book [36J, suggesting a novel approach to better understand the 


traders’ decision making process. One can indeed observe the traders’ immediate reactions 
to price changes in terms of the placement of their orders in the order book. The remarkable 
finding is that the traders’ orders, which may be realised in the future through a trade, 
act similarly to fluid particles ahead or behind a Brownian particle (whose position is the 
mid-price defined as the average of the best bid and ask prices) in the physical fluid analogy. 
Hence the FDT expresses a remarkable relationship between the fluctuations of the mid¬ 
price and an effective drag force exerted by limit orders at different locations in the order 
book. These properties were found to be well described by a Langevin equation. 

These results provide a bridge between the conceptual foundations of the random walk 
theory and its generalisation to both Brownian particles in physics and stock market prices 
in finance. The random walk model of financial prices has been since Bachelier thought to be 
the consequence of the no-arbitrage condition (called the Martingale property in mathemat¬ 
ical finance) and the observable embodiement of the efficient market hypothesis. In contrast, 


our model 


36| in terms of a financial Brownian particle in the layered order book fluid and 


the existence of the fluctuation-dissipation theorem provides a novel picture in which the 
random price walks are, like pollen particles in a fluid, the result of the “collisions” with the 
limit orders that are continuously jiggling, appearing and disappearing. 

Unravelling the physical underpinning of fluctuation phenomena in term of the underly¬ 


ing random wa 
systems [hZ. 38 


k behaviors provides both intuitive and deeper understanding in a variety of 
. Here, we use this novel understanding of financial price dynamics to inves¬ 


tigate the conditions under which the standard mathematical description 


39] of financial 
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price dynamics using Wiener process (the continuous limit of random walks) and stochastic 
Ito calculus hold. Informed by the underlying discrete nature of the limit orders acting as 
effective embedding fluid particles, we introduce a financial Knudsen number, with the pur¬ 


pose of testing the validity of the continuous mathematical models of finance 


40]. In section 


2, we present the data analysis of the foreign exchange data set. In section 3, we introduce 
the Knudsen number of the market. In section 4, we relate the Knudsen number observed in 
the historical time series to the periods when the market price exhibits ballistic-like motions. 
Section 5 concludes. 


II. PRELIMINARY ANALYSIS OF THE DYNAMICS OF LIMIT ORDER BOOKS 

We have analyzed the databse of the Spot FX markets data provided by EBS (ICAP), 
which is known as one of the biggest FX markets mainly for inter-bank dealers. It in¬ 
cludes the order book information of USD/JPY, EUR/USD, and EUR/JPY during the 
week beginning on March 14, 2011. In the data, every order is recorded every millisec¬ 
ond, and a minimum price unit is fixed at 0.001[YEN] for USD/JPY and EUR/JPY, 
0.00001 [USD] for EUR/USD. In the following, the minimum price unit is represented by Ax 
(thus Ax = 0.001 [JPY] for USD/JPY and EUR/JPY, Ax = 0.00001 [USD] for EUR/USD.) 

To analyze the financial time series, we use an event time rather than calendar time. An 
event time is counted after a transaction occurs. After the fc-th transaction, the event time 
t is thus kAt (k = 1,2,...), where At is the average waiting time between two successive 
transactions. 


A. Definition of variables for data analysis 

Let us now present the definitions of the relevant variables. At time t{= kAt), the number 
of limit orders on the price axis x is given by N(x,t)(> 0). The minimum value and unit 
of N(x) is 1 million dollars in USD/JPY, denoted by An (1 million euro in EUR/JPY and 
EUR/USD). 

We represent the best price on the positive (+) side, that is the lowest sell price in the 
order book, as x + (t), and the best price on the negative (—) side, as x~{t) for the highest 
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buy price, and the market price is defined as the mid-price 


x(t) = {x + (t) + x (t)}/2. (1) 

In the time interval [ t — At,t], the velocity of the price (i.e. price change per unit time) is 
defined as 


v(t) = {x(t) — x(t — At)}/At. 


( 2 ) 


The depth of the market is introduced to describe the position of limit orders against the 
best price on its side. At time t, on the j(— +/—) side, the position x of a limit order is 
represented as Qi(t). The depth 7 (Q\t) of the market of a specific limit order with respect 
to the best price at time t — At is defined as 


= sign(j) 


/ Q j (t) - x j (t - At)\ 

V Ax ) ’ 


(3) 


where the depth is counted as multiples of a minimum unit Ax, and sign(+) = 1 and 
sign(—) = —1. In cases when there is no need to show the position x and time t, j^QAt) 
is simply described as 7 , When we need to specify the side j, 7 is shown as 7 T The sign 
in the — side in Eq. ()3l) is needed to make positive the depth of the market. The number 
of limit orders at depth 7 is represented by N^ / ( , y,t). Next, the cumulative volume of limit 
orders from depth 0 to 7 along the market axis depth is defined as 

U(7 = (4) 

7'=0 


Finally, we consider the changes of configuration of limit orders at depth 7 along the 
market axis. In the time interval [t — At, t], the number of limit order changes on the j side 
at position is defined as 


AN^(Q j , t),t) = N(Q j ,t) - N{Q j , t - At) , (5) 

where the 7 dependence of AiV 7 in the l.li.s. is obtained from relationship (|3]) between 7 
and QK Similarly, the change of the cumulative number of limit orders up to 7 is defined 
by 

7 

AE 7 (7(Q^,f),f) = ^A7V T ( 7 ',i) . (6) 

7'=0 
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Fig {H explains the variables defined above. In FigOJa), at time t(= kAt), a configuration 
of the order book is pictured along the horizontal axis x and the vertical axis N(x,t). Blue 
and red circles represent limit orders in the — side and + side. The blue and red dotted 
line represents the limit orders boundary on the — side and + side, that is, x~(t ) and x + (t), 
at time t. In FigJTJb), a configuration of the order book is pictured along the horizontal 
axis 7 and the vertical axis iV 7 ( 7 ,t). Thus, Fig{T^b) is derived from FigOJa) by changing 
the horizontal axis from price x to the depth 7 of the market. The numbers shown below 
the horizontal axis are the values of 7 at each position. In FigOJc), a configuration at 
t + At is shown on the horizontal axis 7 and the vertical axis iV 7 ( 7 , t + At). Black and grey 
circles represent the newly injected limit orders in the time interval [t, t + At], White circles 
represent canceled or executed limit orders. Black and gray circles correspond to 7 > 0 and 
7 < 0. In FigOJd), the change of configuration of the order book between t and t + At is 
shown on the horizontal axis 7 and the vertical axis AiV 7 ( 7 ,t). The changes of limit orders 
between Fig{l](b) and (c), which are shown along the horizontal axis 7 , are encoded in the 
circle colours at each depth 7 . 

The grey circle at 7 = —1 in Fig{l](c),(d) requires some explanation. This negative 7 
corresponds to a price position Q~ = x~(t ) + Ax at t + At falling inside the previous 
interval defined by the best bid and best buy orders at time t. The appearance of this new 
order at time t + At corresponds to the relations N(Q~,t ) = 0, N(Q~,t + At) = 1 and 
AiV 7 ( 7 (Q-,t),t) = 1. While the depth of the new order 7 (Q~,t) is given by 7 (Q~ } t) = 
—{(x _ (t) + Ax) — x~(t)}/Ax = —1 according to definition Eq.Q, it is measured as 7 = 0 
at time t + At. Generally, when 7 (Q\t) < 0, a newly injected limit order is placed lower 
for j = + and higher for j = — than the previous best price in the corresponding ask and 
bid regions. 


B. Properties of the cumulative limit orders as a function of depth 

Let us investigate some properties of the cumulative limit orders up to depth 7, 1 ^( 7 t) 
defined by Eq. (f4]), on both sides j = + and —. Figj2](a) shows the time evolution of 
Fy( 7 + ,t) (where t(= nAt)) for 7 = 0,10 2 ,10 3 (red, green, blue respectively). Although the 
time series for 7 = 0 has many pulse-like peaks, those for 7 = 10 2 and 10 3 seem to behave 
like random walkers. To judge the stationarity of these time series, the power spectra of 


these time series are investigated. In Figj2^b), the power spectra are shown in log-log scale 
(with the same colour code as in Figj2](a)). Approximating a given spectrum by the power 
law S(co) ~ 1 /uj a , recall that a value a larger than 1 in the limit of small uj (large times) 
diagnoses non-stationarity. For 7 = 0, we observe that a is close to 0, and the corresponding 
time series can be regarded as stationary. In contrast, for 7 = 10 2 , and 10 3 , a ~ 2, signalling 
non-stationary time series similar to random walks. This non-stationarity implies that the 
notion of an average limit order structure as a function of depth is meaningless. 

The fact that the number of limit orders at depth 7 = 0 is stationary implies that limit 
orders placed at the best bid and ask prices can be characterised by a well-defined distri¬ 
bution and time-dependence structure. In contrast, the fact that the cumulative numbers 
of limit orders up to large depths follow dynamics similar to random walks, and are thus 
non-stationary, implies to a first approximation that orders are put as random additions, 
cancellations as well as executions without apparent aggregate strategy to keep the total 
order book from accumulating orders as would be the strategy of a market maker trying to 
manage and limit the size of its inventory. 

Although the time series of the cumulative numbers of limit orders (Eq. ([!])) are nonsta¬ 
tionary as shown in Figj 2 ](b) , we checked that the increments defined by eqs. (J5]) and ([ 6 ]) are 
stationary in a weak sense (first-order moment and auto covariance do not vary with time). 
This allows us to study the configuration changes of limit orders caused by newly injected, 
canceled, and executed orders at depth 7 in relation to price changes in a given time interval. 
Let us define the coarse-grained cross correlation function between two variables A(t ) and 
B(t ) by 

Cor(A(t ; kAt), B(t ■ kAt)) = — !—(( A(t ; kAt) - (A))(B{t-, kAt) - ( B ))) , (7) 

O'AO'B 

where (A) and a a represent the average of variable A and the standard deviation (and 
similarly for variable B), and A[t\ kAt ) is given by 

k 1 

A{t; kAt ) = A (t + sAt) . (8) 

where s and k are integer numbers. The coarse graining operation (JHD performed over the 
time scale kAt enables one to remove zig-zag and noisy behaviors at very short time scale to 
extract the robust large-scale correlation coefficients (J7J) . Note that A{t + kAt; kAt ) comes 
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after A(t] kAt) as we do not use the same data points for estimating the two coarse-grained 
variables. 

Substituting v(t;kAt) and AAC^yl, t; kAt) for A(t;kAt) and B(t;kAt) in Eq.((7]), the 
cross correlation coefficients between price changes and the changes of the numbers of limit 
orders at depth 7 are shown in FIG|3ta,b,c) (a: USD/JPY, b: EUR/USD, c:EUR/JPY). 

Blue and red circles corresponds to the — and + side. For these three markets, we observe 
that configuration changes of limit orders have a common functional form of their correlations 
with price changes. In particular, the cross-correlation functions change sign at some critical 
value 7 C defined as follows: 

7 c = In f {7 > 0 : Cor(v(t ), AY 7 (y, t)) < 0, Cor(v(t ), AY 7 (y + 1 , t)) > 0} (9) 

7 “ = Inf{ 7 > 0 : Cor(v(t ), AfV 7 ( 7 , t )) > 0, Cor(v(t ), A 1 V 7 (7 + 1, t )) < 0}. (10) 

FIG13](a,b,c) show that 7 + and 7 “ are essentially undistinguishable, confirming a symmetric 
behavior of the structure of the limit order book in the buy — side and sell + side. 

Substituting v(t\ kAt), A VC,, (7 J , f; kAt) for A(t\ kAt) and B(t\ kAt) inEq.([7|, FIGl3](d,e,f) 
shows the estimated correlation functions between A VC,, (7 J , f; kAt) and v(t\ kAt), where the 
blue and red colour corresponds to the — and + side (d: USD/JPY, e: EUR/USD, f: 

EUR/JPY). One can observe that the maximum of these correlation functions Cor{y{t ; kAt), AVC,( 7 J , t; kAt 
occurs at the previously defined critical values 7 +“ ([9lfT0]h Thus, locating the peaks of 
Cor(v(t] kAt), AV/( 7 U t; kAt)) provides a convenient and robust way for the numerical 
estimation of from the data: 

7 / = sup{ y - 7 > 0 : \Cor{v{t] kAt), AV/(y 7 , f; kAt))\}. ( 11 ) 

The dotted vertical lines in FIGj3])d,e,f) show the location of y c obtained from (fill) : 

(7cA7c + ) = (18,18) for USD/JPY, ( 7 - 7 +) = (17,18) for EUR/USD, ( 7 ^, 7 +) = (22,24) 
for EUR/JPY, where k is set equal to 20. 

These results quantifying the correlations between price increments and changes of the 
number of the limit orders as a function of depth express a collective trend following behavior 
of traders. As price goes up, they tend to requote their limit orders at slightly higher price 
than the best bid. 
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III. KNUDSEN NUMBER OF FINANCIAL MARKETS 


A. Financial Brownian particle in the layered order book fluid 


Previous works have already mentioned an analogy between the order book configuration 
and evolution on the one hand and a fluid of interacting particles on the other hand. In 
particular, Refs. 41 44] note similarities between physics and finance [yj], namely between 
particle collisions in physics and transactions recorded in the order book. In Bak et al.’s 


view 


view 


42], orders are particles and transactions are collisions. In Challet and Stinchcombe’s 
43l] , orders are particles, submission of limit orders are particle depositions, cancellation 


of limit orders represents evaporation, and transactions correspond to particle annihilation. 
Ref. 


has pushed such qualitative analogies to the quantitative level. We introduced 


a description of the dynamics of the order book of financial markets in terms of a Brownian 
particle embedded in a fluid of incoming, exiting and annihilating particles, as illustrated by 
FigSl The Financial Brownian Particle (FBP) is represented by the yellow disk, surrounded 
by the fluid of particles corresponding to the different limit orders existing at time t. The 
vertical dashed line colored with green and orange indicates the position y c Ax from the best 
price in the — and + side. The range x~ — 7 “ Ax < x < x + Ax is called Inner. The range 
x > x + + 7 +Ax, x < x~ — 7 “ Ax is called Outer layer. The intervals x + < x < x + + 7 +Ax and 
x _ — 7 “ Ax < x < x~ are called “interaction ranges” within which the FBP interacts with 
its surrounded particles. The yellow circle delineates the space occupied by the Financial 
Brownian particle, whose size is {x + — x _ }/Ax. 


Figj4](b) shows a configuration change occurring between t and t + At. The newly created 
particles on the side 7 > 0 and 7 < 0) are represented in black and grey. The annihilated 
orders are shown in white. Upward and downward colored arrows indicate the increase and 
decrease of the number of particles in each range, which is supported by the statistical anal¬ 
ysis presented previously in FIGJ3] Figc) shows the new FBP position and configuration 
of the surrounding particles, with in particular the update of the market depth 7 associated 
with the changes of the positions of x + or x~ from t to t + At. 


Within this financial Brownian particle in the layered order book fluid model 36], the 
positive correlation Cor(v(t ; kAt), AlV^y- 7 ', t; kAt )) shown in FIGj3](a,b,c) (a: USD/JPY, b: 
EUR/USD, c:EUR/JPY) for 7 < 7 “ (blue curves on the left side of the black dotted line) is 


11 








interpreted as resulting from a force that pushes the price towards the + side. Symmetrically, 
the negative correlation Cor(v(t) kAt), AiV 7 ( 7 J ", f; kAt)) for 7 < 7 + (red curves on the left 
side of the black dotted line) is interpreted as resulting from a force that pushes the price 
towards the — side. The change of sign of Cor(v(t) kAt), AAi/ 7 / f; kAt)) when 7 crosses 
7 /’“ corresponds to a reversal of the force exerted by the change of limit orders that plays 
the role of a drag resistance against the direction of the price change. 

B. Definition of the mean free path 

Using the above formalism of the FBP, we now relate the variations of the number of 
particles to the motion of the FBP characterised by its velocity v(t). The change in the 
number of particles belonging to the inner layer that is caused by newly created particles is 
noted c*(f)(> 0) and by newly annihilated particles is noted a*(f)(> 0). The change //(f) in 
the number of particles in the inner layer (hence the index i) on the j side is defined as 

7c 

/f (0 = X>Jv 7 (y,o ( 12 ) 

7=0 

= 4 ^)-^)- ( 13 ) 

In FigU c"(f) corresponds to the particles in the green inner layer. c/(f) corresponds to 
the particles in the orange inner layer. These added particles are specifically identified by 
the black colour. As depicted in FigJU a“(f) and a/(f) corresponds to the particles that are 
removed among the green and orange particles in the inner layer. These removed particles 
are represented in white in the inner circle. Then, //(f) is the total flow of particles from f 
to f + At. 

As shown in Fig|3](d,e,f), the cross correlation coefficient between // (t) and v(t) takes 
the largest possible (absolute) value compared to correlations between v(t) and changes in 
particle numbers that would be counted up to values 7 different from 7 C . In fact, the variable 
fi(t) that has the highest cross correlation coefficients with v(t) is derived from //(f), with 
/ = +,-© as 

fi(t) = fr(t)-f+(t). (14) 

Optimizing the coefficients a and (3 of the ansatz af~{t) + /3//~(f) to maximise correlation 
with v(t) gives results that are statistically not significantly better than expression (Tl4]) 
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(which corresponds to a = —/3 = 1). 

Using Eq.(jHJ), we next investigate the relation between the coarse-grained variables 
fi(t\ kAt) and v{t\ kAt). Figj5](a,b,c) shows the existence of a clear linear relation between 
these variables, which can be expressed as 


v(t\ kAt) = L(t ; kAt)fi(t; kAt ) + rj(t; kAt) 


(15) 


where rj(t] kAt) embodies independent noise terms. Here, we have use k = 20. The slope 
L(t] kAt) of the regression has the meaning of a transport coefficient, which is measured in 
the unit [Ax/An). 

Within the physical picture of a Brownian particle embedded in a fluid of incoming, 
exiting and annihilating particles, L can be interpreted as the mean free path of the FBP in 
the order book fluid. Indeed, in molecular dynamics, the mean free path L , collision time r 
and short-time “ballistic” velocity v (between collisions) are related by the formula [45] 


v = 


T 


(16) 


Since there are /, : (t; kAt) collisions in a time interval of duration kAt, as obtained by aver¬ 
aging over a time interval of duration kAt) (|8|) . the typical collision time (i.e. time between 
two collisions) is 


kAt 

fi(t-kAt) 


(17) 


Putting HTTP in (fl6|) yields expression (TT5|l (without the residual term), when counting time 
in units of kAt. 

Estimating L(t]kAt) from the regressions shown in Figj5](a,b,c), we obtain respectively 
L(t ; kAt) = 0.38, 0.34,1.49 for USD/JPY, EUR/USD, EUR/JPY. Thus, the mean free path 
depends on the market. 

As we will show below, financial markets often exhibit asymmetric behavior. This moti¬ 
vates us to generalize (fT5]i into two distinct linear relationships for the + and — sides: 


v^(t; kAt) = V (t; kAt)fl (t; kAt) + kAt) , 


(18) 


where j is either + or —. 
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C. Definition of the Knudsen number Kn 


All financial markets are characterised by discrete price increments Ax and time stamps 
(corresponding to the minimum precision time of 0.001 sec for our forex data). The ques¬ 
tion therefore arises whether the standard financial mathematical formulation in terms of 
continuous stochastic processes [40] hold and under what conditions. A similar question 
is often posed in physics concerning the conditions for the application of the continuous 
Navier-Stokes equation of fluid dynamics, given the fact that the underlying fluid is made 
of discrete molecules. In Physics, this question is addressed by introducing the Knudsen 
number Kn, defined as the ratio of the mean free path of the fluid particles to a character¬ 
istic molecular scale. When Kn is sufficiently smaller than 1, the continuous limit is a good 
approximation of the dynamics. 

Similarly, our model of a financial Brownian particle in the layered order-book fluid 
offers the possibility of defining and estimating a corresponding Knudsen number, which 
characterises a given market. We thus obtain the novel possibility to address quantitatively 
for the first time the question of whether the use of continuous stochastic processes to model 
financial price dynamics is justified. 

We have introduced the mean free path U in expression As a characteristic scale, 
it is natural to use the representative length scale 7 ^ characterising the interaction range, 
shown in Fig. ([4]). Then, the Knudsen number Kn of a given financial market, considering 
possible asymmetric cases, is defined by 


• L j 
Kn J ~ — 

7 c 


A symmetric version of the financial Knudsen number is given by 

Kn 


1 / L + L~ 

2 \7 C + + 7iT 


(19) 


( 20 ) 


For its empirical determination, we estimate L(t\ kAt ) by the least square method in the 
time window [t — (S ■ k — l)Ai, £], which contains S points. In turn, the characteristic scale 
7 c is found to be stable as a function of time, so that we can use a single value for the 
whole period according to the procedure explained from FIGJ2](d,e,f). Fig0a,b,c) showed 
an average value of the mean free path L(t), but this does not mean that L(t) is constant. 
We find in fact that L(t] kAt ) in shorter periods (involving S data points) depends on time, 
so that Kn(t ; kAt) becomes a time-dependent variable. 
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D. Condition for the validity of the continuous time formalism: Kn(t]kAt) < 0Kn 


Figj6]shows the time evolution of the position x(T) of the FBP together with the Knudsen 
number Kn(T ) estimated as described above using k = 4 and S = 100 in three markets 
as a function of calendar time T. We use calendar time to compare markets in which 
transactions do not always occur at the same time. The three markets are USD/JPY, 
EUR/USD and EUR/JPY. The same color code is used for x(t) and Kn(t ) for each currency 
pair. The straight horizontal lines in Figj6](d) corresponds to the average level ( Kn(T )) of 
the Knudsen number over the time period shown. The average values of the Knudsen 
numbers are respectively {0.046,0.039,0.14} for USD/JPY, EUR/USD, EUR/JPY. The 
black dotted horizontal line shows the level 0.1, which is the typical threshold value below 
which the continuum limit is usually considered to be valid. 

A first conclusion is therefore that the continuous limit seems in general adequate for 
the USD/JPY and EUR/USD currency pairs but is more questionable for the EUR/JPY 
pair. For the EUR/JPY pair, this is due to the fact that the position of the FBP (mid-price 
value) fluctuates under the influence of a small number of particles (limit orders), which 
makes discreteness relevant. 

But even for the USD/JPY and EUR/USD pairs, one can observe that there are times 
when Kn(t ) jumps up, signalling periods when the continuous description of the price dy¬ 
namics becomes incorrect. Consider Fig J7] comparing the time series of the USD/JPY ex¬ 
change rate x(t) (position of the FBP) in panel (a) and of its Knudsen number in panel (b). 
The time period covers a regime (depicted in red) in which the price is quite volatile and 
the Knudsen number is above the 0.1 threshold, and another regime in green in which the 
price has low volatility and the Knudsen number is well below the 0.1 threshold. Figj7{c,d) 
magnify the dynamics of x(t) in these two periods. 


E. Relation between mean free path and particle density 

We now empirically demonstrate that high Knudsen numbers occur when the number 
of particles (limit orders) in the interaction range is small, in agreement with the physical 
intuition for the breakdown of the continuous approximation. The total number P(t) of 
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particles in the inner layer on the j side is given by 

7c 

( 21 ) 

7^=0 

We average P(t) over S data points, 

s-k-i 1 

{P(t\ kAt)) s = V /'(« - «'A()fcA«, (22) 

t'=0 

where P is an integer index. We shall also use this time-average over 5 past time instants 
for other variables later. FIGl8](a) demonstrates that the Knudsen number Kn is on average 
inversely proportional to the average number of particles (( I + )s + (I~)s )/2 over the two 
sides averaged over S = 100 data points, 

Kn(t ; kAt) ~ 2k ((/ + (P kAt))s + kAt))s}) 1 , (23) 

for the three currency pairs USD/JPY, EUR/USD and EUR/JPY). In expression ([23]) . we 
use expression (120|) for Kn(t\ kAt) in the symmetric case and thus average the number 
(( I + )s + (I~)s)/ 2 of particles over the two sides. We estimate n = 0.72 by the least 
square method. Given a threshold 6x n (usually set at 0.1) above which the continuous 
approximation is deemed unreliable and discrete effects start to be important, then the 
condition for the validity of the continuous approximation 

Kn{t ; kAt) < 0 Kn (24) 


translates into 


-({/ + (f; kAt)) s + (/ (<; kAt )) s ) > n/0 Kn , 


(25) 


i.e., the average number of particles in the inner layer on each side should be at least 
k/0ku = 7.2 (for Oku = 0.1). Summing over both sides, this implies that at least 14 
particles should be present in the inner layer for the continuous approximation to hold. 


We note that Ref. 


46 ] found a similar relationship for equities between the time-averaged 


number of limit orders at the best quotes and the price changes. Indeed, expression (fl5|) 
together with (1231) implies that v(t; kAt) ~ (( I + (t ; kAt))s + (/“(t; /cAf)^})” 1 . 

Relation (1231) can be generalized to include asymmetric cases as follows: 

c3 


KK (t; kAt) 


K J 


{I j {t ; kAt))s' 


(26) 
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Figj9] shows that Eq. (l26j) holds for sides + and — with k + = 0.71 and k = 0.69, showing 
an almost symmetry behavior. 

Putting together ([T9]) and (I26j) yields 


U (t; kAt) ~ k? 


/ (, Pjt^kAt )) 

V 7c 


s 



(27) 


The variable {P(t] kAt))s/^ c is the ratio of the number of particles in the inner layer to the 
size of that layer. Hence, it is the (linear) density p of particles in the inner layer. Thus, 
equation (1771) writes that the mean free path defined in equation (1151) is inversely proportional 
to the particle density, which is exactly what is expected from the kinetic theory of fluids 


45|. Recall the simple argument to derive this inverse relationship. In a cylinder of length 


L and cross-section a, there are pLcr particles. By definition, the mean free path L is such 
that there is one particle to collide with, if the cross-section of the collision is a. Hence, L is 
determined by pLa ~ 1, which leads to L ~ 1 /pa, which has the same form as (1771) . with the 
identification of p with (P(t; kAt))s/lc- This confirms further the validity of the financial 
Brownian particle model in the fluid of limit orders. Moreover, the coefficient P can be 
interpreted as the inverse cross-section for the collision of the FBP with its surrounding 
limit order particles. 

Putting all this together, Eq. (fT8j) can be generalised to give the dependence of the velocity 
on the j side as 


v J (t; kAt) = lc K ■ fl (t; kAt) + rf (t; kAt) , (28) 

Vi [t] kAt)) s 

where rffakAt) is random noise. Eq. ([28]) stresses that the velocity of the FBP on the j 
side of the inner layer is inversely proportional to the total number of molecules on that side 
in the inner layer. This expression captures the result that, when there is an asymmetric 
change of the number of particles on both sides of the FBP, the total velocity, {u + (t; kAt) + 
v~(t] kAt)}/2 will take an asymmetric value causing a directed motion. In the next section, 
we document examples of such asymmetric market states. 
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IV. ASYMMETRIC PARTICLE DEPLETION RATES, KNUDSEN NUMBERS 
CRITERION AND LARGE VELOCITY OF THE FBP 


A. Asymmetric particle depletion rates 


The generalisation of expression ()28l) (which also makes relation (1T51) more precise) to 
account for asymmetric particle distributions on both + and — sides reads on average 


(v(t; kAt))s = -{(v (t; kAt)) s + (u + (t; kAt)) s } 


7 C • 


■{ft (t-,kAt)) s - 


7c • K 


■{ft(t-,kAt)) s ) , (29) 


2 \{lr(f,kAt))s' Jt (I+(t-kAt)) s 

where 7 / and k? are nearly the same for j — + and j — — . The residual term r ?- 7 is averaged 
out. 

Let us consider the case when the velocity v~ takes a positive value caused by transactions 
or cancelations that totally deplete the number of particles in the inner layer. We want to 
quantify the average speed needed to wipe out the number of particles in the inner layer. 
For this, we measure the ratio of the change (//(£; kAt))s per unit time of the number of 
particles in the inner layer to the total average number (i|(£; kAt))s of particles in the inner 
layer, both on the j side 


\‘(t,kAt)=W t ’ kAt))s 


(30) 


m(t-,kAt))s 

Comparing with expressions (|28li and (l29]h one can see that A J (f;/cAt) is proportional to 
the time-averaged velocity on the j side, up to the constant parameters, 7 / and /v 7 . 

By definition, V{t\kAt) is the rate of change of the number of particles in the inner 
layer on side j. In the continuous limit, this translates into jA = A. For A < 0, it thus 
takes In 2/1A| ticks to halve the number of particles in the inner layer. When the number of 
particles at a certain price level is totally exhausted, the price jumps to the next level. The 
typical time scale for this to occur is ~ 1 /A J '(i; kAt). 

FigJTOl shows the two-dimensional joint probability density distribution of A - and A + . The 
approximate elliptic shape qualifies a bivariable Gaussian distribution, with a negative cross¬ 
correlation —0.58 visible by the tilt of the long axis of the ellipse along the anti-diagonal. 
This means that A~ and A + tend to take opposite signs. When particles are piling up on 
one side, they tend to disappear on the other side. This embodies a transient dominance 
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of buying orders and cancelations. The asymmetric directional motion of the FBP often 
occurs due to an imbalance of the rate of change between the two sides, which is mirrored 
by a corresponding asymmetry of the changes of the Knudsen numbers on both sides. 

Since large asymmetric depletion rates (imbalance of A - and A + ) are often leading to 
strong directed price moves, it is useful to quantify how frequently this occurs. As a bench¬ 
mark, Fig JIUf bl shows the cumulative distribution of 6 defined as the sum over all instances 
in which both A ~ < 9 and A + < 6 occur simultaneously. We find that, in 5% the time series, 
A - and A + are smaller than 0 O .O5 = —0.044. At these rates, the number of particles decays 
by a factor 2 in approximately 6.8 ticks (10.9 seconds on average). We are however mostly 
interested in the asymmetric cases when either A - < 9 or A + < 9 but not both together, 
expressing a strong asymmetric particle depletion on one side. 


B. Joint conditions on asymmetric particle depletion rates and Knudsen numbers 
to detect abnormal market regimes 

Owing to the negative correlation coefficient documented in Fig JIOf a). we expect more 
often the depletion rate to be high on one side while low or average on the other side. 
Fig JTTT a) illustrates this point by showing the locations in red in panel (a) of the USD/JPY 
exchange rate when A + < 6 b.o 5 = —0.044. Fig JTTT bl illustrates this point as well as in cyan 
in panel (b) when A - < 6 0 05 = —0.044 and simultaneously Kn + (resp. Kn~ ) is larger than 
^Kn = 0.1. For this analysis, the above conditions are deemed valid at time t for an analysis 
in the corresponding time interval [t,t — (S ■ k — l)Af], with k = 2 and S = 100 . Fig fTTT c) 
shows the time dependence of the asymmetric Knudsen numbers Kn + and Kn~ with the 
black horizontal dotted line indicating the threshold level 9ku = 0 . 1 . FigJTTTd) presents the 
corresponding time dependence of the rates A + (f) and A ~{t) of particle changes on each side. 

Around t = 125,000, the USD/JPY exchange rate dropped violently. This fall was 
preceded, accompanied and followed by a very large depletion rate on the — side (A_ < 
—0.044) as shown by the presence of the cyan line in Fig JTTlf b) and by a large increase of 
both Knudsen numbers Kn + and Kn~ as shown in Fig JTTT ch This shows that the number 
of particles in the inner layer decreased due to the increase of transactions and cancelations. 
And A~ (in cyan) became negative, and often below the threshold 00.05 = —0.044, while A + 
(in red) remained positive, indicating that a strong imbalance in the depletion speed of the 
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number of the particles in the inner layer occurred in favour of the — side. After t = 127, 000, 
the USD/JPY exchange rate rebounded, with an strong increase in the depletion rate on 
the + side (A + < —0.044) as shown by the presence of the red line in Fig JTTT b) around 
t = 127,000. 

Figffia presents the particle configurations during the abnormal times detected in FigHIl 
Fig fl2lf a) shows the whole time series of the USD/JPY exchange rate x(t) (black) together 
with the flash crash in red. Fig Jl 2 f b) plots the time series x~(t) (blue) and x + (t ) (red) in the 
time interval represented in red in panel (a). Fig fl2lf ci ).(cA and (C 3 ) show the prohles of the 
cumulative number of particles at the best price at the three points indicated in Fig JT 2 f b) 
and located respectively at times t = 124000 (cl) , t = 126000 (c 2 ), and t = 126803 (c 3 ). 
We show respectively U 7 ( 7 \t) in blue for j = — and in red for j = +. Fig JTW dJ shows the 
cumulative number of particles V 1 ( n f~,t) (blue) and U 7 (7 + ,t) (red) as a function of 7 , in the 
limit order book. 

At the time t = 124000 (c 3 ) when the exchange rate begins to drop, the pile of particles 
in front of the FBP exerted a resistance against the downward move. During the period 
when the exchange rate was dropping, and focusing on the specifc time t = 126000 (c 2 ), 
the number of particles on the — side decreased to about 1 / 8 th the number found at time 
t = 124000 (ci), reflecting the downward price motion towards the — side that caused a huge 
annihilation of the particles piled in the — side. Finally, at time t = 126803 (c 3 ), the number 
of particles on the — side totally vanished. In other words, there were no traders on the — 
side. The total number of molecules on the — and + side shown in Fig Jl2f dl as a function 
of the depth in the limit order book indicates that the profile of the number of particles on 
the — side in log-scale remained approximately the same during the duration of the flash 
crash and its rebound. This confirms that the changes in the number of particles strongly 
contribute to the price formation during the crash. Theoretically, the Knudsen number on 
the — side reaches infinity at the bottom t = 126803 (c 3 ). This is a very rare situation when 
there is no buyer of dollar in the market and no one can sell dollar any more. This halt of 
the market function occurred in the middle of a working hour. The above analysis shows 
that high values of the Knudsen numbers and the development of asymmetric properties of 
the change of particle numbers on each side of the inner layer can detect this peculiar event 
at its early stage. 
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V. CONCLUSION 


Starting from the description of the dynamics of the order book of financial markets 
in terms of a Brownian particle embedded in a fluid of incoming, exiting and annihilating 
particles [36], we have generalised this model by presenting a model of the velocity on each 
side (buy and sell) independently. The improved model builds on the time-averaged number 
(Jj) of particles in the inner layer and its change f t per unit time (the index i stands for “inner 
layer”). The inner layer has been defined by specific robust properties of the layering of the 
particles as a function of the order book depth, revealed by the correlations between price 
velocity and change in the number of particles. This allowed us to introduce the Knudsen 
number of the financial Brownian particle (FBP) motion and its asymmetric version (on the 
buy and sell sides). Not being considered previously, the asymmetric Knudsen numbers are 
crucial in finance in order to detect asymmetric price changes. We even found rare regimes 
when no particle stands on one side, so that the corresponding asymmetric Knudsen number 
is in principle infinite, as shown in FigJl2l Such situation is associated with an extraordinary 
market regime, such as a flash crash, at which the market stops functioning. 

By measuring the Knudsen numbers, we have shown that it is possible to characterise 
the conditions for the market dynamics to be correctly described by a continuous stochastic 
process. Even though the continuum formulation is usually applied without questions for 
large liquid markets such as the USD/JPY and EUR/USD exchange rates, we have shown 
that there are regimes when the Knudsen numbers are so high that discrete particle ef¬ 
fects dominate, such as during market stresses and crashes. Moreover, we found that the 
EUR/JPY market operates most of the time at rather large Knudsen numbers, so that the 
continuous formulation should be amended. 

We also found that the Knudsen number is inversely proportional to the number of 
particles (limit orders) in the inner layer. This confirms the condition for the breakdown 
of the continuous formulation, i.e. when the number of particles is too small for a local 
average behavior to be defined. Indeed, small numbers of particles translate in large Knudsen 
numbers, i.e., the mean free path is larger than the typical particle sizes. If observable, 
the number of limit orders in the inner layer can be a good indicator of the stability and 
directional movement of financial markets. 

Finally, we documented an imbalance of the depletion rates of the number of particles on 
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the buy and sell sides, with the clear occurrence of high Knudsen numbers associated with 
violent directional price changes. This indicator is more useful than the popular symmetric 
volatility measure. With these metrics, we could detect the direction of the price motion at 
the early stage while volatility is blind to the price direction. 
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FIG. 1. Schematic representation and definition of variables encoding the information contained 
in the order book, (a) A configuration of the order book at time t is shown on the horizontal price 
axis x and the vertical axis N(x,t). Blue circles correspond to limit orders on the — side and red 
circles correspond to limit orders on the + side. The dotted line coloured in blue indicates the 
best price on the — side and red indicates the best price on the + side, (b) The configuration in 
(a) is represented on the horizontal 7 3 axis and vertical t) axis. Under the horizontal axis, 

the numbers indicate the depth values, (c) A configuration of the order book at time t + At is 
shown on the horizontal and vertical axis AT( 7 , t + At). Newly injected orders are coloured in 
black for 7 > 0 or grey for 7 < 0, while canceled or executed orders are white, (d) The difference 
between the number of limit orders at t + At and t is shown on the horizontal 7 J axis and vertical 
AT 7 ( 7 , t) axis. 
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FIG. 2. Study of the stationarity of the number of limit orders as a function of time: (a) 
cumulative number Vy(j,t) of limit orders up to depth 7 as a function of time t for 7 = 0 (red), 
7 = 10 2 (green) and 7 = 10 3 (blue), (b) Power spectra S^(uj) of the three times series of panel (a) 
are shown as a function of angular frequency w, with the same colour code as in panel (a). At long 
time scales (w < 10 —2 ), we can observe that S 7 =o(uj) becomes flat, i.e. with an exponent (defined 
in the text) converging to 0 , diagnosing a stationary behavior. In contrast, both S' 7=10 2(w) and 
5 7=10 3 (u;) have their exponent a « 2, which is larger than 1, thus expressing a non-stationarity 
behavior of (7 = 10 2 , t) and (7 = 10 3 ,i). 
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FIG. 3 . Study of the relationship between the dynamics of the changes with time of the number 
of limit orders and the price changes: (a) Cross correlation function Cor(v(t), AIV 7 (7, £; kAt)) for 
USD/JPY as a function of 7, showing the critical depth 7 c defined in the text with the black 
dotted line; (d) Cross correlation function Cor(v(t), AV/(7, f; kAt)) for USD/JPY as a function of 
7, where the critical depth 7 C is indicated with the black dotted line; Panels (b,e) are the same as 
panels (a, b) for the EUR/USD exchange rate. Panels (c,f) are the same as panels (a, b) for the 
EUR/JPY exchange rate. Blue and red circles corresponds to the — and + side. 
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FIG. 4 . Financial Brownian Particle (FBP) represented by the yellow disk with size {x + — x~ + 
Ax}/Ax., surrounded by the fluid of particles corresponding to the different limit orders existing 
at time t. (a) A configuration of the FBP and surrounding particles is represented along the 
horizontal position x of the order book price and the vertical axis N(x,t ) counting the particles 
(orders). The vertical dashed line colored with green indicates the position 77 Ax from the best 
price in the — side. Orange one indicates the position 7 +Ax from the best price in the + side. 
The black dotted line shows the position x(t) of the centre of the FBP. (b) A configuration change 
N(x,t + At) occurring between t and t + At. N(x, t + At) gives the new numbers of limit orders as 
a function of price x. The newly created particles on the side 7 > 0 and 7 < 0 ) are represented in 
black and grey. The annihilated orders are shown in white. Upward and downward colored arrows 
indicate the increase and decrease of the number of particles in each range, (c) New FBP position 
and configuration N(x, t + At) of the surrounding particles at time t + At Dotted coloured vertical 
line have moved from their previous positions at time t (panel (a)) because of the modification of 
the best prices, x~(t) and x + {t). Correspondingly, the interaction range is also updated, which 
revise the colours of the particles. 
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FIG. 5. Empirical evidence for Eg. (11511 : (a) For USD/JPY, scatter plot of u(t; kAt ) as a function 
of fi(t]kAt), with k = 20. The median is the thick blue dot and the error bars are the inter- 
quartiles. The fitted dotted line is estimated by the least square method, (b) Same as (a) for 
EUR/USD. (c) Same as (a) for EUR/JPY. 
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FIG. 6. Time evolution of the position x(T) of the FBP together with the Knudsen number 
Kn(T), using k = 4 and S = 100 in three markets as a function of calendar time T. (a) USD/JPY, 
x(T) (red); (b) EUR/USD, x{T) (green); (c) EUR/JPY, x(T) (blue); (d) Kn(T) for the three 
currency pairs with the same colour code as in (a),(b), and (c). The colored dotted lines show the 
average values of Kn. The black dotted horizontal line shows the level 0.1. 
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FIG. 7. Typical time evolution of the USD/JPY exchange rate and of the corresponding Knudsen 
number: (a) Exchange rate time series x{t) (position of the FBP); the red is enlarged in panel 
(c) and the blue is enlarged in panel (d). (b) Estimated Knudsen number corresponding to the 
time series shown in panel (a), (c) Enlarged time series of the exchange rate x(t), corresponding 
to an estimated Knudsen number above 0.1. (d) Enlarged time series of x(t), corresponding to an 
estimated Knudsen number well below 0.1. 
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FIG. 8. Dependence of the Knudsen number as a function of the number of particles in the inner 
layer, (a) Average Knudsen number Kn(t ; kAt ) in the symmetric case estimated using expression 
( 1201 ) as a function of the average number (( I + )s + (I~)s)/2 of particles over the two sides, with 
S = 100. We use the same color code as in Figj6]to refer to the three pairs of currencies: USD/JPY 
(red), EUR/USD (green), EUR/JPY (blue), (bi) Median and quartiles of the data points in (a) 
together with the best fit with expression (flOli (dotted blue line) k = 0.72. (b 2 ) log-log plot of (bi) 
illustrating the inverse relationship (|23D . 
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FIG. 9. Dependence of the Knudsen number as a function of the number of particles in the inner 
layer estimated separately for the + and — sides: (ai) Median and quartiles of Kn~ (t; kAt) as a 
function of ( I~(t ; kAt))s (with S = 100) for the same three currency pairs as in Fig. [HI The best 
fit with expression (1261) (dotted line) gives k~ = 0.69. ( & 2 ) log-log scale of (ai) illustrating the 
inverse relationship (1261) . (bi) same as (ai) for the + side, with k + = 0.71. (b 2 ) log-log scale of 
(bi). 
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FIG. 10. (a) Joint histogram of the rates of change of particle numbers A - and A + in the inner 

layer exhibiting a negative cross-correlation of —0.58. (b) Cumulative distribution of 0 defined as 
the normalised sum over all instances in which both A ~ < 6 and A + < 0 occur simultaneously. The 
horizontal dotted line shows the 5% confidence level, which is reached at 6 = —0.044. 
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FIG. 11. Time dynamics of Knudsen numbers and rates of the change of particle numbers in the 
inner layer around a flash crash for the USD/JPY exchange rate on March 16, 2011 (GMT+0). 
On the horizontal axis, time is in multiple units of tick time, (a) Locations in red of the USD/JPY 
exchange rate when A + < #o .05 = —0.044 and simultaneously Kn + is larger than O^n = 0.1. (b) 
Locations in cyan of the USD/JPY exchange rate when A - < #o .05 = —0.044 and simultaneously 
Kn~ is larger than Qku = 0.1. For given values at time t, the evaluation of these variables is 
performed in the corresponding time interval [t,t — (S ■ k — l)Af], with k = 2 and S = 100. (c) 
Time dependence of the asymmetric Knudsen numbers Kn + (t ; kAt) (red) and Kn~(t ; kAt) (blue) 
with the black horizontal dotted line indicating the threshold level &Kn = 0.1. (d) Corresponding 
time dependence of the rates A + (f) (red) and A ~(t) (blue) of particle changes on each side with 
the black horizontal dotted line indicating the threshold level $o .05 = —0.044. 
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FIG. 12. Particle configurations during the abnormal times detected in Figllll (a) Whole time 
series of the USD/JPY exchange rate x(t) (black) together with the flash crash in red. (b) Time 
series x~(t ) (blue) and x + {t) (red) in the time interval represented in red in panel (a). (ci),(c 2 ) 
and (C 3 ) show the profiles of the cumulative number of particles at the best price at the three 
points indicated in Fig flYlf bl and located respectively at times t = 124000 (ci) , t = 126000 ( 02 ), 
and t = 126803 (C 3 ). ^(y - ,^) is shown in blue and V y ('y + ,t) is shown in red. (d) Cumulative 
number of particles Fy( 7 _ ,t) (blue) and ^(7 + ,t) (red) in logarithmic scale as a function of 7 in 
logarithmic scale in the limit order book. 
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